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Design criteria to minimize modal dispersion have been found for 
a broad class of practical, multimode, circular-symmetric, isotropic, 
optical fibers having any numerical aperture and any profile dispersion 
(which is a function of the derivative of the index with respect to the 
wavelength). The impulse-response width of these fibers, the rms width 
of the impulse response, the optimum profiles to minimize those widths, 
and the sensitivity to profile departures from ideal are found to be 
surprisingly simple closed- form generalizations of previous results that 
are mostly applicable to fibers with small numerical aperture and 
constant profile dispersion. The minimum impulse-response width of 
the optimized fiber is a function only of its numerical aperture and 
consequently is independent of the index profile and of the profile 
dispersion. 



I. INTRODUCTION 

Circular-symmetric, multimode, optical fibers intended for large 
communication capacity must have low modal dispersion and this is 
achievable by the quasi-complete equalization of the group velocities 
of all modes 1 (or rays). This equalization depends critically both on the 
refractive-index profile and on the profile dispersion of the fiber. The 
profile dispersion is defined in Section II, but here it is enough to know 
that it is related to the derivative of the index with respect to the wave- 
length. 

To understand better the objectives of this paper, let us first review 
some recent evolution of thoughts linking the index profile and the 
profile dispersion of a fiber to the pulse broadening caused by modal 
dispersion. 
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Gloge and Marcatili 2 showed that if the numerical aperture (NA) of 
the fiber is arbitrary but the profile dispersion is negligible, there is a 
family of fibers— for which the dielectric constant profiles decrease ra- 
dially according to power laws — that is important for two reasons. The 
first reason is that the family encompasses a wide variety of easy-to-make 
fibers (step-index, quasi-parabolic, etc.) possessing the unique property 
that the group velocity of each mode is a function only of its propagation 
constant; this drastically simplifies the analysis. The second and more 
important point is that for an almost parabolic power law of the dielectric 
profile, a fiber with small NA has the very narrow impulse response 
needed for high-speed communication. 

Olshansky and Keck 3 extended these results in a very important way 
by showing that if the profile dispersion is constant across the core, 
narrow impulse response is achievable in small NA fibers by a simple 
modification of the exponent of the dielectric-constant profile's power 
law. 

In many cases, though, the two requirements — smallness of NA and 
constancy of profile dispersion— are not satisfied. For example, to in- 
crease the coupling efficiency to incoherent sources and to decrease 
microbending losses, 4 fibers with large NAS 8 are being made. They are 
heavily doped and, particularly if the doping element is boron, the profile 
dispersion may not be a constant 6 - 7 as a function of the radius. Similar 
lack of constancy may occur in fibers that are doped with several mate- 
rials for the purpose of improving optical or mechanical properties. 8 
Arnaud and Flemming 9 - 10 have calculated the impulse response for these 
fibers, treating the variable portions of the profile dispersion as a per- 
turbation. Using a numerical method Arnaud 11 has also calculated the 
pulse spreading in a multimode planar fiber with arbitrary index profile 
and profile dispersion. 

In this paper, we extend the previous results by finding, within the 
WKB approximation, a surprisingly simple closed-form description of 
the modal dispersion in a broad class of circular-symmetric isotropic 
fibers which have arbitrarily large NA and arbitrary profile disper- 
sion. 

The gist of our paper is in Sections II and III. In Section II, the profiles 
of the fibers belonging to the group are defined and their impulse-re- 
sponse widths are calculated. In Section III, the optimum index profile 
required to minimize the impulse-response width is determined together 
with the sensitivity of this response to departures of the index from 
optimum. The rms of the impulse response is the subject of Section IV. 
In Section V, some approximate results about the influence of index 
profiles on the minimization of the rms width of the impulse response 
are derived and conclusions are drawn in Section VI. 
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II. FIBER PROFILE AND WIDTH OF ITS IMPULSE RESPONSE 

We start by looking for the dielectric-constant profile of a circular 
symmetric isotropic fiber such that, as in Ref. 2, the group velocity of 
each mode is only a function of its propagation constant. In the process, 
we will also find the width of the impulse response of that fiber. 

The initial point is a WKB approximation 2 that relates the propagation 
constant 0(X) of a mode characterized by the radial and azimuthal wave 
numbers n and v to the free-space propagation constant k = 2ir/\, the 
refractive index n(r, X) of the fiber and the radial coordinate r, via the 
integral 

,-i rA (i) 

TT Jr\ r 

where 



p(r, X) = V(k 2 n 2 -0 2 )r*-»* (2) 

and ri and r> are two neighboring turning points that make the radical 
zero and between which most of the field of the mode is concentrated. 
It is useful to redefine 

n 2 = n 2 (l-F) (3) 

2 = k 2 ni(l-B), (4) 

where n \ is the index on axis and the profile function F(r, X) is zero on 
axis, is an arbitrary function of r and X within the core {r ^ a), and is 
2A(X), a function only of X in the cladding (r ^ a). Similarly, the mode 
parameter B varies between zero for the lowest-order mode and 2A(X) 
for the modes whose phase velocities coincide with that of a plane wave 
in the cladding.* With these definitions, the radical (2) becomes 

p = V(kn l r) 2 (B-F)-u 2 . (5) 

The group velocity of a mode (or ray) is introduced by taking the de- 
rivative of both sides of (1) with respect to the free-space wavelength 
X, 



rN'-ftiS-'^-a];*-* 



where 



X dn\' 



/ A an\\ 



(6) 



(7) 



* Similar but not identical profile function and mode parameter have been introduced 
previously in the literature. 12 
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is the group index on axis and 

p{rA) = ^FTx {8) 

is a generalized version of the profile dispersion parameter introduced 

in Ref. 12. 

The derivative dB/d\ in (6) can be expressed in terms of the group 
delay t of the mode by taking square roots and derivatives on both sides 
of (4). The result is 



dkm T ' 2N x VT=Bd\ 

in which T, the flight time on axis, that is, the delay of a plane wave in 
a medium of group index N\ and length L, is related to the velocity of 
the light in free space c via 

r-£&- do) 

C 

The substitution ofdB/dX from (9) into (6) yields the integral 

r[ 1 -* =B ?-'( 1 -9];*- a (11) 

This integral was solved in Refs. (2) and (3) by assuming p constant and 
F = 2A(r/a)", a power law, with a an arbitrary constant. To lift these 
restrictions and still solve (11) exactly, the following self-evident ex- 
pression is introduced: 



C' 2 ^-dr = p(r 2 ,\)-p(ruX) = 0. 
Jr\ or 



(12) 



This integral becomes less obvious and very useful when the derivative 
dp/dr is performed with the help of (5), yielding 

r-\ B _ F _rJlldr = (13) 

J n L 2dr] p 

Combining (13) with (11), we arrive at a general expression 

'-^? .0-5)7* 



B r r-> , r dF\ Fr , 

■ — ) — dr 



(14) 



J rx \ 2FdrI p 
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that like (11), relates the group delay t of a mode characterized by its 
mode parameter B (or propagation constant ft), and its azimuthal mode 
number /' (hidden in p) to the profile function F and profile dispersion 
p. This expression is valid for any circular-symmetric fiber with isotropic 
dielectric and, in general, still cannot be solved exactly. However, if a 
particular family of fibers is considered that satisfies the condition 

r dF 

l+ Y F Tr 

=D(\), (15) 

2 

D being an arbitrary function of X, the seemingly formidable right-hand 
side of (14) is reduced to l/D and the group delay of the mode charac- 
terized by the mode parameter B is 

t=T vf^- • ,16) 

These last two equations are the basic results of the paper. Equation 
(16) says that t, the group delay of a mode (or ray) is only a function of 
the mode parameter B and the dispersion parameter D. More important, 
the group delay is independent of the mode number (which means that 
modes with the same propagation constant have the same delay), it is 
independent of the profile function F, and it is independent of the profile 
dispersion p. Equation (16) is used in the following sections to study the 
impulse response of the fiber. 

On the other hand, eq. (15) is the recipe for the design of the fiber 
whose time response is given by (16). It can be solved in several ways 
depending on the control that the fiber designer has over F and its de- 
rivatives with respect to X. The least demands on these functions occur 
if the fiber is designed to operate at one wavelength only. Then, in (15), 
D becomes a constant, p is only a function of r, and the partial derivative 
of F is reduced to an ordinary one. Without introducing new symbols 
for D, F, and p, the simplified design formula is 

r dF 

1 + VFTr 

=D. (17) 

2 

This equation in turn can be solved in two ways. One way consists in 
prescribing the profile function F to satisfy, perhaps, requirements 
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different from modal dispersion. Then, the profile dispersion p must 
be tailored to satisfy (17). 

For example, assume the dielectric profile depicted in Fig. la. The 
profile function is 

2a/— )"' for 0<r<a o 
/ =. f W (18) 

-\ for a <r<a, 



where 



and the inequalities 



/ai\«2/(«-«) 



a\ < a 

OC\> CX2> 



(19) 



(20) 



guarantee that the profile looks indeed like that in Fig. la. 

Substituting (18) in (17) and assuming for the dispersion parameter 
D a value D that optimizes in some sense the impulse response of the 
fiber, the required profile dispersion turns out to be 



P = 



2 _2±^1 for r<fl0 

°° (21) 

2 _2±^2 for ao<r<a 

Dn 



and is shown in Fig. lb. 

This is an interesting example not only because it clearly illustrates 
the power of the theory even to design optimized fibers in which the 
profile function and dispersion are discontinuous, but also because it 
may be of practical interest. For example, by using an index-increasing 
dopant for r <a and an index-decreasing dopant for r > a the NA of 
the fiber can be increased, keeping at the same time its optimum modal 
dispersion. 

In the other way of solving (17), the profile dispersion p as a function 
F is assumed to be known, from experiment, and the index profile must 
be found from the integration of (17) that yields 

J»2A dp 
7 ^ TT^' (22) 

f [2-D(2-p)]F 

This result will be used in a more general way later, but if for the time 
being we prescribe p to be a constant P , the profile function results: 
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where 



(b) 
Fig. 1— (a) Dielectric profile (solid line), (b) Profile dispersion. 

F=2aQ'\ (23) 

a = D(2 - P ) - 2. (24) 



This last equation establishes the relation between the dispersion pa- 
rameter D of the fiber introduced in this paper and the a value so widely 
used in the literature 2 -' 5 for fibers with constant profile dispersion Po. 
It follows from (17) that only if p is a constant, is the profile function F 
a power law (23). 

MODAL DISPERSION IN OPTICAL FIBERS 55 






The two solutions described require only the control over the profile 
function F and its first derivative with respect to X. But suppose that 
the fiber designer has also control over the second derivative. Then, to 
increase the range of wavelengths over which the fiber operates with low 
modal dispersion, he could demand, for example, that at a wavelength 
not only 

D = D Q 
but also, as proposed by Kaminow and Presby, 13 

^=0. (25) 

dX 

This requires the simultaneous satisfaction of (17) and 

d /r dF\ 

^l*±=- Do (26) 

dp 

d\ 

derived from (15) and (25). It is this last equation that implies the control 
over the second derivative of F with respect to X. 

It can easily be extrapolated that control over higher derivatives 
permits even further demands on D. In fact, if all the higher derivatives 
were controllable, D(X) could be chosen arbitrarily and the profile F 
would be the solution of the linear partial differential equation of first 
order (15) subject to the conditions of being zero at r = and 2A(X) at 
r = a; the result is well known 14 from a mathematical point of view, but 
of limited importance from a practical point of view. 

III. MINIMIZATION OF THE IMPULSE-RESPONSE WIDTH AND ITS 
SENSITIVITY TO ERRORS IN THE PROFILE 

The impulse-response width is determined from (16) by finding the 
difference between flight times of the slowest and the fastest modes (or 
rays) for any given value of the dispersion parameter D. It is simple to 
find that the minimum time spread, T min , between those modes occurs 
if D is chosen 

D =l + V1-2A. (27) 

In fact, the modes characterized by B = and B = 2 A are the slowest and 
arrive at the end of the fiber at 

U. ■ T, (28) 

while the modes characterized by B - 1 - Vl - 2A are the fastest and 
arrive at 
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2(1 - 2A)'/« . 

'""" T i + VT^IK ( ' 

The time spread between them is the minimum impulse-response 
width 

[1-(1-2A)*/']' . 

r m in - t m ax ~ *min - i j + ^ _ gA WW 

Therefore, fibers with the same A have the same minimum impulse- 
response width r m in, independently of their index profiles and profile 
dispersions, provided that they satisfy the design equation (17) with D 
substituted by the optimum value D (27). 
If 

A«l, 

(27) and (30) become 

Z) = 2 - A (31) 

rmin s V T = °- 61 A2 » s/km (32) 

8 

for iVi- 1.46. 

To find the sensitivity of the impulse-response width to departures 
of the index profile from optimum, we calculate the ratio rhmm between 
the response width t for 

D = (1 + 8)D , (33) 

where 

5«1 

and the minimum response width r min occurring for D = D . After some 
straightforward calculations, 



7", 

and for 



1 + ^ (1 + \/T^2A) 2 [1 + (1 - 2A) 1 ' 4 ] 2 
8A 



(34) 



-^-/l + W. (35) 

Tmin V A / 

It is known that the impulse-response width is indeed very sensitive to 
the choice of profile and more so for smaller A. If 5, the fractional de- 
parture of D from its optimum value, is equal to A, then the pulse width 
is nine times larger than r mm . 

The main results in the last two sections have been extended by Ar- 
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naud to optimize modal dispersion in fibers with noncircularly symmetric 
profiles. 15 

IV. THE RMS OF THE IMPULSE RESPONSE 

From the point of view of the maximum information-carrying capacity 
of a fiber, more significant than the impulse response width is its rms 
width a, 16 since l/4a is the repetition rate at which pulses can be trans- 
mitted with a reasonable loss penalty at the receiving end. 1 

Let us calculate first the impulse response W(t) and then its rms width 
a assuming that: 

(£) The energy of the infinitely narrow impulse fed at one end of the 

fiber is equipartitioned among all modes. 
(») All modes attenuate equally. 
(ill) The number of modes is so large that the discrete pulses arriving 

at the receiving end can be replaced by a continuum. 

The impulse response, then, is the rate of change of the number of 
modes reaching the end of the fiber, 

W{t) =4~ (" mtt \(",t)dv, (36) 

at Jo 

and its rms width is, by definition, 



1 J Jo 



1/2 

W{ti)W{t 2 ){ti - t 2 ) 2 dtidt 2 y 



SSo" mt 



(37) 



)W(t 2 )dt l dt 2 



To calculate W(t), the value of n given in (1) is substituted in (36) and 
the integration along v is carried through yielding 

Wit)= ^l± C r " (B -F)rdr. (38) 

4 dt Jo 

The integral measures the energy arriving at the end of the fiber as a 
function of time and, since each contribution must be positive, the largest 
value that F can reach is B. Therefore, the upper limit r B is the value of 
r that makes 

F(r, X) = B(X). (39) 

The explicit value of r B depends on the choice of fiber design. If the 
profile function F is prescribed, then r B is obtained by solving (39). If, 
on the other hand, the profile dispersion is prescribed, then 

J»2A dp 

B [2-D(2-p)]F m 
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follows from (22). 

Now the derivative with respect to time in (38) is carried out. The 
derivative of the integral is equal to the integral of the derivative since 
the terms that should contain the integrand times the derivatives of the 
limits are zero. Consequently, 

*(,)•&£&£ (41) 

The reader interested in the explicit impulse response must substitute 
B in this expression with its time-dependent value obtained from 
(16). 

Replacing W(t) in (37) and also substituting the explicit value of t 
from (16), the rms width of the impulse response results in 

,M lA Iv^ VT^ | | (42) 



\ ) ) r'irjdxdy 



where x and y are dummy variables and r x and r y are given by (40) once 
B is substituted either by x or by y . It is easy to recognize in (42) that if 
A « 1 and D ^ 2, the parenthesis is of the order of A 2 and a is propor- 
tional to A 2 T. 

Unlike the simple impulse-response width, the rms width <r and the 
optimum value of D that minimizes it are dependent on the profile dis- 
persion p and the profile function F. In general, the exact value of a and 
its minimizations must be found numerically, but we push the analysis 
a little further in the next section where some simplifying assumptions 
are made. 

V. APPROXIMATE RESULTS FOR RMS WIDTH OF THE IMPULSE 

RESPONSE, ITS MINIMIZATION, AND ITS SENSITIVITY TO PROFILE 
ERRORS 

Within the family of fibers described in the previous sections there 
is a large group of particular importance that encompasses many of the 
available fibers today. This group has small NA and its profile dispersion 
is almost constant with respect to r. To introduce these properties, we 
assume 

A « 1, (43) 

then the profile dispersion is expanded in power series of the profile 
function F, 

I F \ « 

P = £ Ps (—) , (44) 

s =o \2A/ 
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and since F is a function of r, the near invariance of p with r implies 

-£y « i. ( 45 > 



i,m^) 



5.1 Prof lie function 

Carrying this simplifying assumption to (22) by keeping only first 
powers of P s (s > 0), the profile function results: 

where 

a = D(2 - P ) - 2 (47) 



and D is still an arbitrary number. 

If the profile dispersion is constant, the summation in (46) disappears; 
then, and only then, will the profile function follow a pure power law. 

5.2 Minimization of the rms width of the Impulse response and Its sensitivity 
to profile errors 

We want to find <x min , the minimum rms width of the impulse response 
possible, and Dj, the optimum dispersion parameter for which <r mm is 
achieved. The optimum profile is obtained by substituting D with Di 
in (46). We are interested also in finding the sensitivity of o to small er- 
rors in the profile. 

To achieve these purposes o 2 is expanded in a power series about Di, 
and only the first three terms are kept, 

da 2 (D-Di) 2 dV 

The derivatives are to be taken at D - D\. Since by definition a 2 passes 
through a minimum of D = D\, the equation 

— = at D = D } (49) 

dD 

serves to determine the optimum dispersion parameter D\. 

From (42), (48), and (49), we obtain with the help of (43) and (45) 
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^Rirl^H 



(50) 



(1 + H) 1/2 
a min = T(Atf) 2 a + sH){l + 4H){1 + 5H) U2 (51) 

-2- = A/1 + ^ " Dl * + — ) 2 (1 + 3g>(l + 5H)> (52) 
a min V VD.AH 1 + 2H/ 

where 

- (s 2 + s + 6)ff 2 + 8#+2 

2 " 2 s ?i J {[(a + 2)H + l] 2 - ff 2 }{[(« + S)H + l] 2 - H 2 } 

and 

// = 1 - P . (54) 

The optimum value of the dispersion parameter D\ is close to 2. The 
profile function that maximizes the information-carrying capacity of 
the fiber is obtained by substituting D with Di in (46). The dispersion - 
profile terms of order higher than zero appear in (50) only in the sum- 
mation 2 and are multiplied by A. Therefore, their contribution is small 
indeed and is neglected in (51) and (52). It is kept in (50) because, as will 
be seen later, small errors in the profile affect substantially the value of 
a. For 2 = 0, the substitution of (50) in (47) yields the same optimum 
aofRef. 3. 

Consider now a m i n » the minimum rms width of the impulse response. 
From (51), we might be tempted to conclude that/f = 1 — Pq should be 
made small to decrease a mm . However, the number of modes of the guide 
derived from (38) with the help of (46) is 

" = ff) 2 TT7T 

Therefore, if the number of modes of the fiber is to be kept constant, ff m i n 
can be decreased by making 

{l + H) 5 ' 2 



(i + 3//)(i + 4H)(i + my' 2 



small and this is achieved by choosing H large, not small. 

The following table contains the minimum rms per unit length of fiber, 
0min/^» and the concomitant pulse repetition rate L/Ao m \ n for different 
values of// as derived from (58), assuming N\ = 1.46. 
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H a m JL PRR = TT- 

Ms/km Mb -km/s 



1 0.14 A 2 1-79/A 2 

2 0.16A 2 1.56/ A 2 
0.18 A 2 1.38/ A 2 

For H > 1, the pulse repetition rate is fairly insensitive to the value of 
H. ForH = 1 and A = 0.01, the pulse repetition rate is ~18 Gb/km/s. This 
information-carrying capacity is only 33 percent smaller than that of the 
"ideal profile" reported by Cook. 17 

Let us turn now to the sensitivity of the rms width to errors in profile 
(52). Again, for H > 1, this result is insensitive to the value of H; in- 
deed 






(56) 



For H = 1 and (D - Di)/D lt the fractional departure of D with respect 
to the optimum D x equal to A, the rms width a is about 8.4 times wider 
than a min . As in the case of the pulse width, the rms width is very sensi- 
tive to profile errors. 

A fiber designed to minimize the rms width (D = Di) has only 30 
percent more information-carrying capacity than a fiber with the same 
A designed to minimize the impulse-response width (D = D ). 



VI. CONCLUSIONS 

For a vast class of circular-symmetric fibers made of isotropic di- 
electrics, simple and fundamental design criteria that minimize the 
impulse-response width due to modal dispersion at one wavelength have 
been found. This minimum width (30) is only a function of the NA and 
the time of flight along the axis. Therefore, if properly designed, a fiber 
with arbitrary profile dispersion has the same minimum impulse-re- 
sponse width as another fiber with the same NA and no profile dispersion. 
Their information -carrying capacity is about 1.4/A 2 Mb • km/s. The fiber 
engineer has a substantial freedom of choice to reach that optimum 
design: the profile dispersion may be arbitrarily chosen but then the 
index profile is uniquely determined by (22); or symmetrically, the index 
profile may be arbitrarily chosen and then the profile dispersion must 
satisfy (17). Only if the profile dispersion is a constant does the optimum 
dielectric profile that minimizes the impulse response follow a power 
law. 
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The profile dispersion entails the first derivative of the index with 
respect to the wavelength. If the second derivative can be controlled, then 
the minimization of the impulse-response width can be achieved at two 
neighboring wavelengths. This broadbanding of the fiber response can 
be expanded even further if higher derivatives are under control. 

The width of the impulse response is very sensitive to errors in the 
fiber design. A fractional error of A in the dispersion parameter of the 
fiber makes the response about nine times wider than the minimum as 
seen from (35). 

Only a marginal increase of about 25 percent in the information- 
carrying capacity of the fiber is achieved if the rms width of the impulse 
response is minimized instead of minimizing the impulse-response 
width. 
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